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^ '. 1 Introduction 

\Q • We are closing in on neutron stars both observationally and theoretically. Observa- 

^ . tionally, a number of masses (M) and a few radii (R) have been measured as well 

^ ! as a number of other properties. Theoretically, modern equation of states (EOS) are 

^ I more reliable due to precision measurements of nucleon-nucleon interactions, detailed 
calculations of binding energies of light nuclei and nuclear matter which constrain 

' three-body forces, inclusion of relativistic effects, improved many-body and Monte 

^ . Carlo methods, etc. 

C\J ! Ultimately, we can exploit the one-to-one correspondence between the EOS and 

O ■ 

o ■ 

' Observing a range of neutron star M and R thus reveals the EOS (e.g., pressure 

• P versus density p) of dense and cold hadronic matter. Possible phase transitions 

^ ! from nuclear matter to quark matter (either can also undergo superfluid transitions 

I at certain densities and temperatures), hyperon matter, kaon or pion condensates, 

' etc., would also be revealed by an anomalous/kinky function P{p) and M{R). The 

higher the order of the transition is, the smoother will M{R) be and very accurate 
observations will thus be necessary. On the other hand, just one accurately measured 
neutron star mass and radius would already constrain the EOS significantly. Infor- 
mation on the EOS at high baryon densities and the presence or absence of phase 
transitions could guide us in solving QCD after decades of unsuccessful attempts. 

In the following I shall give a brief account of the present status on neutron star 
observations and theory referring to ^, ^, Q for longer reviews. Subsequently, I shall 
attempt to recount the most likely possible phase transition in dense nuclear matter 
with emphasis on quark matter and its possible color superconducting states, as this 
is most relevant at this conference. Finally, I shall point to important developments 
expected in the near future. 



the mass-radius relation of cold stellar object: 

P{p) ^ M{R) (1) 
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Figure 1: Neutron star masses vs. radius for modern |^ and FPS ||^ EOS and strange 
stars 0. The hatched areas represent the neutron star radii and masses allowed for 
orbital QPO frequencies 1060 Hz of 4U 1820-30 (vertical lines, and for burster 

oscillations of 4U 1636-53 assuming M > IAMq (horizontal lines, |l^) area. Models 
for glitches in the Vela pulsar constrain masses and radii []rT| below the full line. The 
radii of RX J1856-3754 from Refs. 



12, 13 assumes M = 1.37Mq. 



2 Observed neutron star masses 

Only a few masses have been determined from the more than thousand neutron stars, 
that have been discovered so far: 

Binary pulsars: Six double neutron star binaries are known so far, and all of 
them have masses in the surprisingly narrow range 1.36±0.08Mq WM. Neutron stars 



are estimated to have a binding energy of ~ 10% of their mass. Thus ~ 1.5 Mq of 
nuclei are needed to obtain a 1.35 Mq star. It is suspicious that the Chandrasekhar 
mass (maximum mass before gravitational collapse sets in) for the iron core inside 
a large burning star is just around ~ 1.5 M©. It is therefore a tempting conclusion 
that the iron cores are the progenitors of neutron stars and that all neutron stars are 
simply produced with M ~ 1.35 Mq. Similarly, white dwarfs are formed in a narrow 
mass range around M ~ 0.6 Mq whereas their Chandrasekhar mass is M ~ 1.35 Mq. 
The latter mass is probably reached by accretion and is responsible for supernova type 
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SN-la used as standard candles in cosmology. Calculations of supernova explosions 
do, however, indicate that neutron stars should be formed with a wide range of 
masses and so the narrow binary pulsar mass range could be a selection effect in 
forming a double neutron star system. Whereas selection effects often are important 
in astrophysics the contrary is the case in particle physics, where one does not believe 
in accidents but invoke some underlying symmetry. 

Millisecond binary pulsars with white dwarf companions have less precise mass 
determination. Promising progress is reported, e.g., for PSR J0437-4715 |]l5l where a 
pulsar mass of M = (1.58 ± O.18)M0 (error bars are 2a) is found. 

Vela X-1 and Cygnus X-2 are X-ray binary pulsar/burster with high/low 
mass companions respectively. From X-ray pulse delays, optical radial velocities and 
constraints from X-ray eclipse, their masses have been determined. For Vela X-1: 
M = 1.87^;]:?? Mq [g, and for Cygnus X-2: M = 1.8 ± 0.4Mq fl^. 



QPO's are neutron stars emitting X-ray's at frequencies of the orbiting accreting 
matter. Such quasi-periodic oscillations (QPO) have been found in 12 binaries of 
neutron stars with low mass companions. If the QPO originate from the innermost 
stable orbit ||^, Q of the accreting matter, their observed values imply that the accret- 
ing neutron star has a mass ~ 2.3Mq in the case of 4U 1820-30. If not, the QPO's 
constrain the EOS as shown in Fig. 1. 



3 Observed neutron star radii 

The small size of neutron stars makes it very difficult to observe them directly and 
measure their radius. Estimates have been obtained using quite different methods, 
which has the benefit that the systematic errors are also different. 

RX J1856. 5-3754 is our nearest known neutron star. It is non-pulsating and 
almost thermally radiating. It has been studied recently with the Hubble space 



telescope by Walter et al. |T^. Its surface temperature is T ~ 57 eV and its distance 
is, from parallax measurements and circumstantial evidence, about d ~ 61 pc. From 
the measured flux 

F = asBT'R'/d' (2) 



one obtains a radius of Rao = R/yl — 2GM/R ~ 7 km which is incompatible with 
almost any EOS. Kaplan et al. |jl9l have reanalysed the HST data and find only 
half the parallax and thus twice the distance and radius Roo — 15 km corresponding 
to i? ~ 12 km for M = 1.4M0. Its age would be almost a million years which 
is compatible with standard modified URCA cooling. The spectrum is, however, 
suppressed in the optical part as compared to the X-rays. Recent more detailed 
analyses of the spectrum |l^, ^ attempting to model the neutron star atmosphere 



and its absorption, as well as including magnetic fields does not improve the spectral 
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fit. A much better description is obtained from a two temperature model, i.e., a small 
hot spot and a colder area with a larger radius ~ 9((i/61pc) km. 

Gravitational lensing by our nearest neutron star may be observed in June 2003, 
when J1856. 5-3754 will pass within ~ 0.3 arcseconds of a background star as it flies 



through space with proper motion of 0.332 arcseconds per year. According to [21 



(see, however, [0) it should just be possible to measure that the 26.5 magnitude 
background star moves by 0.6 milliarcseconds due to the gravitational field of the 
neutron star. This accuracy requires that the Hubble space telescope is extended 
with the Advanced Camera for Surveys. If possible, its mass could be derived and 
our nearest known neutron star would then also be the first one with both mass and 
radius determined. 

X-ray bursts are thermonuclear explosions of accreted matter on the surface of 
neutron stars. After accumulating hydrogen on the surface for hours, pressure and 
temperature become sufficient to trigger a runaway thermonuclear explosion seen as 
an X-ray burst that lasts a few seconds E^. Assuming that the burst spectra are 



black-body one can from the resulting temperature and measured flux estimate the 
neutron star radius if its distance is known. Often the radius is underestimated 
because only a hot spot emits or the spectrum contains a hard tail. Some bursts do, 
however, give radii of order ~ 12 km with a period of almost constant (Eddington) 
luminosities. 

Quescent neutron stars are non-accreting X-ray binaries of which some are 
emitting thermally. Recently, a few quescent neutron stars have been discovered in 
galactic globular clusters where the distance is known relatively accurately. CXOU 
132619.7-472910.8 in NGC 5139 has T = 66 ± 5 eV and = 14.3 ± 2.1 km (90% 
confidence limit) and similar radii although with larger uncertainties are obtained 
from half a dozen other quescent neutron stars . 

Burst oscillations in the X-ray fiux during the first seconds of the bursts have 
recently been exploited to determine the compactness or redshift M/R. Whereas the 
average amplitude increases with the growing size of the hot spot, the amplitude is 
strongly modulated by the rotational period of the neutron star. The fiux does not 
completely disappear when the hot spot is on the back side due to bending of the 



light in the gravitational field. This way M/R < 0.16 can be extracted |T0| for 4U 
1636-53 and a radius i? > 12 — 13 km is obtained for a M = 1.4M0 neutron star (see 
Fig. 1). Corrections from aberration, doppler shifts, etc. are being investigated. Yet, 
the oscillation analyses is another new promising method by which we can obtain 
neutron star masses and radii. 

Absorption lines in the neutron star photospheres should be detectable with the 
spectrographs on board Chandra and XMM. The gravitational redshift and pressure 
broadening of absorption lines determine M/R and M/R? respectively and would 
thus complement other mass and radius information. First results were, however, 
disappointing. Only in RX J1856. 5-3754 were there indications for one or two lines 
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| T3| but higher resolution is required and will come with future observation time. 



4 Modern equation of states of dense matter 

The relation between M{R) and P{p) is often presented (in particular by observers) 
by a wide variety of curves calculated from "old" EOS thus implying that one can get 
anything from theory. This is very unfair as many old EOS's are inconsistent with 
known nucleon-nucleon interactions and/or saturation density and binding energy of 
nuclei/nuclear matter and should therefore not be taken seriously. 

Modern microscopic EOS's are actually converging ^. The NN interaction is 
now well determined and constrain potential models leaving only minor differences. 
Relativistic effects have been included and current scattering experiments at inter- 
mediate energies will determine the NN interactions at higher momenta relevant for 
higher densities. Three-body interactions can be constrained in order to fit nuclear 
binding and saturation density as well as binding energies of light nuclei up to A < 8. 
Many-body and Monte Carlo techniques are now much more accurate. The resulting 
"modern" EOS are quite reliable up to a few times nuclear saturation densities. 
Above they are expected to break down but can be constrained by causality condi- 
tions. It is boldly predicted that neutron stars in the mass range M ~ 0.8 — l.SM© 
all have radius just around R ~ 11.5 km (see Fig. 1). 

The uncertainties in the EOS at densities p^Spo affect M{R) for the heavy neutron 
stars only. By making the EOS stiffer at high densities in a smooth way, the maximum 
mass can increase up to M^2.0 — 2.3Mq ^ but not much higher due to the causality 
(cs < c) condition. Rotation can increase it by ~ 10% for the millisecond pulsars. 

Phase transitions generally soften the EOS and lower the maximum mass. If 
the M = 2.25Mq mass in 4U 1820-30 stands this would rule out any major phase 
transition and allow only the stiffest EOS. 

Many other phase transitions in neutron stars have been considered. It is expected 
that superfiuid phases of neutrons and protons exist at least in certain density regions 
at low temperatures although have not been calculated reliably for strongly inter- 
acting and correlated nuclear matter. At typical neutron star densities neutrons and 
protons are superfiuid as well due to ^Sq and, in the case of protons, also pairing 
0]. These superfiuid and superconducting components will have drastically different 
transport properties than normal Fermi liquids. Generally the resistance, specific 
heat, viscosities, cooling, etc. are suppressed by factors of order ~ exp(— Aj/T), 
where Aj is the gap of quarks, nucleons or electrons. A superfiuid neutron gas in the 
inner crust is assumed in the description of glitches |Tl| and provides constraints on 
the EOS (see Fig. 1). 

Eventually at very high densities nucleon degrees of freedom must be replaced by 
quark ones but it is not known whether core densities of neutron stars are sufficient. 
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Such quark stars will be discussed in the following section. 

More speculative phases are 7r°, vr" and K~ condensates as well as hyperons 
A,...). In p a condensate of virtual ttO is found in a narrow density interval due to 
strong tensor correlations. The K~ energy can be calculated at low densities and 
a naive extrapolation would lead to condensate at high densities p^4po- However, 
correlations in nuclear matter invalidates such an extrapolation and makes a K~ 
condensate unlikely [0] . Hyperons are found to appear at rather low densities p^2 — 
3po in a number of models |2^. Due to limited information on hyperon-nucleon two- 
and three-body interactions one cannot exclude the presence of hyperons in cores 
of neutron stars but their effect on the binding energy and thus the EOS is minor 
whereas their effect on /ig could be substantial 0. 



Figure 2: Cross section of a ~ IAMq neutron star. The ~ 1 km thick crust consist of 
neutron rich nuclei in a lattice and a uniform background of electrons and, in the inner 
crust, also a neutron gas. The interior of the neutron star contains a nuclear liquid 
of mainly neutrons and ~ 10% protons at densities above nuclear matter density no 
increasing towards the center. Here pressures and densities may be sufficiently high 
that the dense cold strongly interacting matter undergoes phase transitions to, e.g., 
quark or hyperon matter or pion or kaon condensates appear. Typical sizes of the 
nuclear and quark matter structures are ~ 10~^^ m but have been scaled up to be 
seen. 
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5 Quark stars 



Natures marvelous variety of EOS's results in the existence of several different stable 
stellar objects all around one solar mass. These can be ordinary stars, white dwarfs, 
neutron stars, black holes - and possibly also quark stars. Quark stars come in several 
categories depending on the details of the nuclear to quark matter phase transition: 

• The pure quark stars also called "strange" stars consist of up, down and strange 
quarks with electrons to fulfil charge neutrality. Possibly a crust of nuclei is 
suspended above the surface of the quark star. In simple bag model EOS a 
rather low bag constant and strange quark mass is required to make strange 
stars and strangelets. If the SAX J1808. 4-3658 really has i? ~ 6 km based 
on accretion in magnetic fields [0 or RX J1856-1754 has i? ~ 7 km from the 



one-temperature fit in Ref. [|T2[ that would indicate strange stars rather than 
normal neutron stars. 

Hybrid stars have a core of quark matter and a mantle of nuclear matter. The 
quark core size depends on the EOS and vanishes for large bag constants leading 
to a normal (nucleons only) neutron star. 

Mixed stars have a mixed phase of nuclear and quark matter over a range of 
density or radius. The mixed phase appear in two-component systems, where 



the two components: neutrons and protons or up and down quarks [^]. It 
is, however, required that the interface tension is sufficiently small so that the 
surface and Coulomb energies of the associated structures are small ||2^. If not, 
a hybrid star results. 

A mixed phase of quark and nuclear matter has lower energy per baryon at a 
wide range of densities if the Coulomb and surface energies associated with the 



structures are sufficiently small [^, The mixed phase will then consist of two 
coexisting phases of nuclear and quark matter in droplet, rod- or plate- like structures 
(see Fig. 2) in a continuous background of electrons much like the mixed phase of 
nuclear matter and a neutron gas in the inner crust of neutron stars [^. Another 
requirement for a mixed phase is that the length scales of such structures must be 
shorter than typical screening lengths of electrons, protons and quarks. 

In the mixed phase the nuclear and quark matter will be positively and negatively 
charged respectively. /3-equilibrium determines the chemical potentials and densities 
in the coexisting phases. Total charge neutrality 

2 11 

ne = {I- f)np + fi-Uu - -Ud - -ris) , (3) 

where Up is the proton density, determines the "filling fraction" /, i.e. the fraction 
of the volume filled by quark matter. For pure nuclear matter, / = 0, the nuclear 
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symmetry energy can force the electron chemical potential above ~ 100 MeV at a 
few times normal nuclear matter densities. With increasing filling fraction, however, 
negative charged droplets of quark matter replace some of the electrons and fie de- 
creases. With increasing density and filling fraction it drops to its minimum value 
given fie = ml/Ajj, corresponding to pure quark matter, / = 1. 
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Figure 3: Sketch of the QCD phase diagram, temperature vs. baryon density in 
neutron star matter, i.e. charge neutral and in /5-equilibrium containing electrons. 
Hatched areas indicate mixed phases of hadronic matter (HM) and quark matter 
(QM/QGP) as well as the nuclear liquid-gas. Dash-dotted lines indicate melting 
temperatures of the lattices in the mixed phase. Dashed lines separate CSC phases 
that may appear (see text). The trajectory of neutron star core densities during 
formation is shown by dotted line and densities below it exist inside neutron stars. 



6 Color superconductivity 

If quark matter exists in neutron stars or is produced in heavy-ion collisions, a con- 
densate of quark Cooper pairs may appear at low temperatures characterized by a 
BCS gap A usually referred to as color superconductivity (CSC) The appearance 
of a gap through color-flavor locking (CFL) requires the gap to exceed the difference 
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between the u,d,s quark Fermi momenta, which is not the case for sufficiently large 
strange quark masses or for an appreciable electron chemical potential, fig, which is 
present in neutron star matter as discussed in [p2 |. 



In neutron star matter /5-equilibrium relates the quark and electron chemical po- 
tentials 

fJ'd = f^s = fJ'u + fJ'e ■ (4) 

Temperatures are normally much smaller than typical Fermi energies in neutron stars. 

The strange quark masses and electron chemical potentials stress the system in 
the direction of splitting the quark chemical potentials. The pairing interaction prefer 
overlapping quark Fermi surfaces. We shall investigate this competition in detail 
below. The case where interactions between quarks are strong and the effective pairing 
gaps is larger than A^v^|/^e/2 — m^/8/i| was shown in pO| to favor the CFL phase. 



If interactions are weak, the chemical potentials are then related to Fermi momenta 



by fii = ymf +pf. If the strange quark mass is much smaller than the quark 
chemical potentials, Eq. (Q) implies a difference between the quark Fermi momenta 

Pu-Pd = /ie , (5) 

2 

m 

Pu-Ps ^ T, /^e, (6) 

Pd-Ps ^ — , 7 
2/i 

where /i is an average quark chemical potential. Strange quark masses are estimated 
from low energy QCD ~ 150 — 200 MeV and typical quark chemical potentials 
are typically /i ~ 400 — 600 MeV for quark matter in neutron stars 0]. Consequently, 
m2/2/i~ 10-25 MeV. 

The BCS gap equation has previously been solved for pure u,d and u,d,s quark 
matter ignoring electrons and /5-equilibrium and the conditions for condensates of 
dicolor pairs (2CS) and CFL respectively were obtained. There are three 2CS condi- 
tions 

A>|pi-pj|, i^j = u,d,s. (8) 

In bulk quark matter total charge neutrality of quarks and electrons and /^-equilibrium 
require that /ig — m^/4/i. In a mixed phase of quark and nuclear matter the electron 
chemical potential is a decreasing function from the value in pure /3-equilibrium nu- 
clear matter fie ~ 100 — 200 MeV down to that for bulk quark matter (if the cores of 
neutron stars are very dense) /ie — m^/4/i ^5 — 10 MeV. 

We can now give a qualitative picture of the various phases that strongly interact- 
ing matter undergoes as the density increases towards the center of a cold neutron star 
(see Fig. 3). The inner crust will undergo several transitions as the nuclear matter 
and neutron gas mixed phase change dimensionality via nuclei, rods, plates, tubes. 
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bubbles to nuclear matter in which the neutron and protons may be superfluid. If 
quark matter appears at higher densities in a mixed phase with nuclear matter those 
structures repeat again. 

Inside the quark matter part of this mixed phase the CSC phase also changes with 
density or /Xg depending on the size of the effective pairing interaction A between two 
quarks i,j (see, e.g., p3[). 



If the pairing is strong A^v^|/ie/2 — m^/8yu| the CFL phase is favored and no 
electrons appear as was shown in . This condition is not likely to be fulfilled 
in the mixed phase with low quark filling fraction, where fie may be of order 
100 - 200 MeV or larger. 

If /ie/2^A^mg/2/i a number of CSC phases appear. In the beginning of the 
mixed phase (/ ~ 0) only fulfils (||) and we have a 2CS of d,s quarks. A 2CS 
of u,s quarks may, however, compete at larger filling. At the end of the mixed 
phase (/ ~ 1) all (||,|,^ fulfil (||) resulting in CFL. In between a crystalline 
LOFF phase and a CFL phase with a vr" condensate (analogous to the CFL-i^" 
phase in [^) may appear. 

If the pairing interaction is weak or the strange quark mass large such that 
A^m^/4/i a number of different CSC phases such as a 2CS of u,s quarks or 



CFL-K may appear as has been discussed in p^, ^ . 



In the mixed phase gaps may be affected by the finite size of the quark matter 
structures as is the case for nuclei. Also surface and Coulomb energies generally 
disfavors mixed phases. 

The finite temperature extension of the competing phases, calculations of the 
critical temperatures and densities, order of the transitions, etc. for neutron star 
matter should be investigated further. Probably, the superfluid phases undergo a 
second order phase transition to the normal phase with increasing temperature at 
constant baryon density. However, transitions between competing phases might also 
occur as indicated in Fig. 3. 

Temperatures in neutron stars, T^IO^K ~ 10^^ MeV after cooling, are typically 
much lower than lattice melting temperatures, Tmeit — Z'^e^/VJ^a where a is the 
lattice spacing (see Fig. 3). Thus the quark matter structures would be solid frozen 
and the cores of neutron stars would be crystalline and possibly also CSC. Lattice 
vibration will couple electrons at the Fermi surface with opposite momenta and spins 
via phonons and lead to a "standard" BCS gap for electrons. The isotopic masses are 
similar but as densities and Debye frequencies are larger, we can expect considerably 
larger BCS gaps for electrons. The electrical super conduct ivy affect magnetic fields 
through the Meissner effect and magnetic field decay. 
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7 Outlook 



A number of promising developments have been mentioned above that may provide 
new information on neutron star masses, radii, EOS and possible phase transitions 
in high density matter. 

On the observational front more and better masses and radii are determined by a 
number of different methods (thus with different systematic errors). On the theoret- 
ical front the uncertainties in the EOS are reduced by improved two- and three-body 
forces, relativistic effects and manybody calculations leading to socalled modern EOS. 
A number of phase transitions are still possible at high densities and can lead to mar- 
velous structured phases and condensates. 

We can hope for additional information from other directions as well: 
Gamma ray bursters. The discovery of afterglow in Gamma Ray Bursters 
(GRB) allows determination of their very high redshifts {z > 1). They imply that 
GRB have an enormous energy output ~ 10^^ ergs which requires some central engine 



more powerful than ordinary supernovae These could be a special class of type 



Ic supernova (hypernovae) where cores collapse to black holes, or binary neutron stars 
merging, or some major phase transition to, e.g., quark matter ^5]. Also soft GRB 
may be explained by accretion on strange stars . 

Neutrinos from the formation of a "proto-neutron star" will be detected from 
the predicted 1-3 supernovae explosions in our and neighboring galaxy per century. 
In the case of the recent 1987A in LMC 19 neutrinos were detected on earth and 
with the upgraded neutrino detectors many thousand neutrinos are expected. The 
neutrinos can test the SN models, the neutron star EOS and early cooling. During the 
first second rapid cooling takes place and a phase transitions to, e.g., quark matter 
or CSC may occur. This may result in a delayed neutrino blimps |^ . 



Gravitational wave ground based interferometic detectors currently under com- 
mission will improve sensitivity significantly. Detectable candidate sources are inspi- 
ralling binary neutron stars or black holes merging which may be responsible for 
gamma ray bursts. R-mode instabilities in rapidly rotating neutron stars may also 
be detectable |8|, ggf. 

Relativistic heavy-ion collisions are probing the high temperature and low 
baryon density part of the QCD phase diagram at the opposite end to neutron stars. 
The phase transition from hadronic matter to a quark-gluon plasma is searched for 
though the transition may well be a smooth cross over. Some "anomalous" effects, 
i.e. deviations from predictions based on hadronic theory, have been found and are 
currently studied intensively at RHIC. It has been speculated that one might find a 
critical point at high temperatures which would then indicate that the smooth cross 
over changes to a first order transition at higher baryon densities. That would be 
valuable information for neutron stars although it would not tell at which density the 
transition would occur at low temperatures. 
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Superfluidity in cold fermionic atomic systems will, due to the rapid progress 
in cooling and trapping techniques, soon be discovered. The great advantage of 
atomic traps is that one can vary the number of particles, their density and interac- 
tion strength as well as the number of spin states, masses, etc. Measuring the size of 
the pairing gap and varying these parameters provides a testing ground for analytical 
calculations in the dilute or weakly interacting limit This would be very useful 



for gap calculations in neutron and nuclear matter as well as quark matter. 

The future of neutron star observations looks bright as new windows are about to 
open. A new fleet of X- and Gamma-ray satellites have and will be launched. With 
upgraded ground based observatories and detectors for neutrinos and gravitational 
waves |^9[ our knowledge of neutron star properties will be greatly improved. Heavy- 
ion physics at RHIC and LHC may add further to our understanding of the QCD 
phase diagram. 

Thanks to R. Ouyed and F. Sannino for organizing this conference and comments 
on the manuscript. 



References 

[1] J.M. Lattimer and M. Prakash, Phys. Rep. 333, 121 (2000). 
[2] S. Balberg, S.L. Shapiro, |astro-ph/ 0004311 . 

[3] H. Heiselberg and V.R. Pandharipande, Ann. Rev. Nucl. Part. sci. 50, 481 (2000). 

[4] H. Heiselberg and M. Hjorth- Jensen, Phys. Rep. 328/5-6, 237-327 (2000). 

[5] A. Akmal, V. R. Pandharipande, D. G. Ravenhall, Phys. Rev. C 58, 1804 (1998). 

[6] B. Friedman and V.R. Pandharipande, Nucl. Phys. A361, 502 (1981). 

[7] X.-D. LI, 1. Bombaci, M. Dey, J. Dey and E.P.J, van den Heuvel, Phys. Rev. 
Lett. 83, 3776 (1999); 1. Bombaci, |astro-ph/020"T369| . 

[8] W. Zhang, T. E. Strohmayer and J. H. Swank, Ap. J. 482, L167 (1997). 

[9] M. C. Miller, F. K. Lamb and D. Psahis, Ap. J. 508, 791 (1998). 

[10] N.R. Nath, T.E. Strohmayer and J.H. Swank, |astro-ph/0102421| |astra 



ph/0102421| . T.S. Olson, Phys. Rev. C 63, 015802 (2001); |astro-ph/020T099 



CM. Miller, priv. comm. 
[11] B. Link, R. I. Epstein, J. M. Lattimer, Phys. Rev. Lett. 83, 3362 (1999). 



12 



[12 
[13 
[14 

[is: 

[16 

[17; 

[18, 

[19 
[20 
[21 
[22 
[23 



[24 

[25 
[26 
[27; 
[28 
[29 

[30 
[31 



J. A. Pons et al., ApJ 564, (2002), in press. 

M.H. van Kerkwijk, |astro-ph/0110336i 

S.E. Thorsett and D. Chakrabarty, Ap. J. 512, 288 (1999). 

W. van Straten et al., |astro-ph/0108254| . 



O. Barziv et al., |astro-ph/0108237| . M. H. van Kerkwijk, |astro-ph/0001077| 
(2000). 

J.A. Orosz and E. Kuulkers, Mon. Not. R. Astron. Soc, 305, 1320 (1999). 

P.M. Walter and L.D. Mettliews, Nature 389, 358 (1997); P.M. Walter, S.J. 
Wolk, and R. Neuhauser, Nature 379, 233 (1996). 



D.L. Kaplan, M.H. van Kerkwijk, J. Anderson, |astro-ph / 1 1 1 1 7^ . 
V. Burwitz et al., |astro-ph/0109374 . 



B. Paczynski, |astro-ph/0107443| . 

L. Bildsten and T. Strohmayer, Physics Today, Peb., p. 40, (1999). 

E. P. Brown, L. Bildsten and R. E. Rutledge, Ap. J., 504, L95 (1998); |astra 



ph/0 1054051 . 



V.R. Pandharipande et al, Phys. Rev. Lett. 75, 4567 (1995); Phys. Rev. C 63, 
017603 (2001). 

N. Glendenning and J. Schaffner, Phys. Rev. Lett. 81, 4564 (1998). 

P. Weber, |astro-ph/0112058| . N.K. Glendenning, Phys. Rev. D 46, 1274 (1992). 



H. Heiselberg, C.J. Pethick, E.P. Staubo, Phys. Rev. Letts. 70, 1355 (1993). 

CP. Lorenz, D.G. Ravenhall and C.J. Pethick, Phys. Rev. Lett. 70, 379 (1993). 

See, e.g., M. Alford, K. Rajagopal, S. Reddy, P. Wilczek, Phys.Rev. D 64, 074017 
(2001), and refs. therein. M. Alford, |hep-ph/0102"047i |hep-ph/01ldT50| . D.T. Son, 
|hep-ph/ 01082601 . 



K. Rajagopal, P. Wilczek, Phys. Rev. Lett. 86, 3492 (2001). 

K. Splittorff, D.T.Son, M.A. Stephanov, Phys.Rev. D64, 016003 (2001). J.B. 
Kogut, D. Toublan, Phys.Rev. D64, 034007 (2001). 



[32] H. Heiselberg, |hep-ph/99124T9 



13 



[33] J.E. Horvath et al., |astro-ph /0112159| . D. Blaschke, S. Fredriksson, A.M. Oztas, 
lastro-ph/011 15871 . 

[34] T. Piran, lastro-ph/Oliml . 

[35] D.K. Hong, S.D.H. Hsu, F. Sannino, Phys.Lett. B516 (2001) 362-366. R. Ouyed 
and F. Sannino, |astro-ph/0103022| ; |hep-ph/01120T3| . 



[36] V.V. Usov, |astro-ph/01lTli2 



[37] G. Carter, |hep-ph/0111353 



[38] J. Madsen, |hep-ph/01114T7 . 

[39] See, e.g., K. Belczynski et al., |astro-ph/01 1 1452] , for a discussion of LIGO can- 
didates. 

[40] H. Heiselberg et al., Phys. Rev. Letters, 85, 2418 (2000); Phys. Rev. A63 043606 
(2001) 



14 



